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We address the issue whether quasifree single-nucleon knockout measurements carry sufficient 
information about the nuclear interior. To this end, we present comparisons of the reaction prob- 
ability densities for A(e,ep) and A(p,2p) in quasifree kinematics for the target nuclei ^He, ^^C, 
^^Fe, and ^°^Pb. We adopt a comprehensive framework based on the impulse approximation and on 
a relativized extension of Glauber multiple-scattering reaction theory in which the medium effects 
related to short-range correlations (SRC) are implemented. It is demonstrated that SRC weaken 
the effect of attenuation. For light target nuclei, both the quasifree (p, 2p) and (e, ep) can probe 
average densities of the same order as nuclear saturation density po. For heavy nuclei like ^°^Pb, 
the probed average densities are smaller than O.lpo and the (e, ep) reaction is far more efficient in 
probing the bulk regions than (p, 2p) . 
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I. INTRODUCTION 

Single-nucleon knockout reactions from nuclei in 
quasifree kinematics continue to be major source of in- 
formation of the mean-field properties of nuclei. In a 
quasifree A(e, e'p) reaction, a bound nucleon in the tar- 
get nucleus A is subjected to a electron- nucleon interac- 
tion and ejected, thereby leaving the residual nucleus in 
a low- lying hole state. The shape of the measured dif- 
ferential cross sections and the knowledge about the en- 
ergy of the residual nucleus, allows one to determine the 
quantum numbers of the nucleon that was struck by the 
virtual photon. Electroinduced single-proton knockout 
experiments with stable nuclear targets have systemati- 
cally pointed towards the validity but also the limitations 
of the mean-field picture for understanding nuclei jHE]. 

The measured quasifree A{e^ e'p) differential cross sec- 
tions provided evidence for the mean-field picture in that 
the extracted momentum distributions for the bound nu- 
cleons could be modeled with mean-field single-particle 
wave functions. It should be stressed that the measured 
momentum distributions cannot be directly related to 
the single-particle wave functions in momentum space, 
but are distorted in the sense that they are affected by 
the attenuation effects of the nuclear medium on the 
ejected proton. The systematic observation that the ex- 
tracted spectroscopic factors from A{e^ e'p) studies are 
substantially smaller than the predicted single-particle 
level occupancies provide evidence for the limitations 
of the mean-field picture [Sj. The extracted spectro- 
scopic factors are defined as the overall normalization fac- 
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tor between the measured and the computed differential 
cross sections. The computed A{e^ e'p) observables are 
based on models that make assumptions for the electron- 
nucleus interaction, the nuclear wave functions, and the 
effect of nuclear attenuation on the ejected proton. The 
latter effect is often referred to as final-state interactions 
(FSI). The credibility of a reaction model for A{e^ e'p) de- 
pends on its ability to describe for example the differen- 
tial cross sections and polarization observables for a num- 
ber of target nuclei. An interesting question is whether 
the extracted normalization factors can be related to the 
shell occupancy in a model-independent fashion [4, 5 . 

An alternative method to gain access to the mean-field 
properties of nuclei is the A{p^2p) reaction [Q. For the 
study of stable nuclei one may prefer the A{e^ e'p) reac- 
tion as it leads to some reduced parameter dependence 
of the extracted results, like spectroscopic factors. In the 
first place, this is due to the electromagnetic character of 
the interaction vertex. Second, for obvious reasons the 
sensitivity to the modeled nuclear attenuation is larger 
for (p, 2p) than for (e^e'p). Electron scattering experi- 
ments from unstable nuclei are a real technical challenge 
and could for example be performed at an eA collider 
like the one which is on the drawing table for the ELISe 
(ELectron-Ion Scattering in a Storage Ring) project at 
FAIR [7^. In inverse kinematics (i.e. the p{A^2p)A — 1 
process) the A{p, 2p) reaction offers great perspectives 
to investigate the mean-field properties of unstable nu- 
clei [8 . One of the fundamental questions, for exam- 
ple, we need to ask ourselves is how the single-particle 
properties of nuclei evolve as a function of the proton- 
to- neutron ratio. Recent studies [9 with the heavy-ion- 
induced nucleon-knockout reactions ^Be{A^ A — 1)X^ sug- 
gest that the deduced spectroscopic factors for single- 
proton and single-neutron knockout are dramatically de- 
pendent on the asymmetry of the proton and neutron 
Fermi surface. The ^Be{A^ A — 1)X process is extremely 
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surface dominated and the above-mentioned results with 
regard to spectroscopic factors await confirmation with 
a reaction probe which is more efficient in probing the 
nuclear interior. 

The scientific potential of the (p, 2p) reaction in in- 
verse kinematics to study the mass number dependence 
of mean-field properties, is nicely illustrated in Ref. jTOj . 
There, results of (p, 2p) measurements on the eight car- 
bon isotopes 9-1^ C are presented. The measurements 
provide empirical information about the mass-number 
dependence of the weakly-bound and inner-shell protons. 
It is shown that the systematics of the momentum dis- 
tributions, separation energies, and spectroscopic factors 
as a function of the mass number can be studied. The 
(p, 2p) measurements of Ref. [10 are performed at ener- 
gies of 250 MeV/A. In the foreseeable future, quasifree 
(j9, 2p) measurements for protons of several hundreds of 
MeV/A will become feasible with high-energy heavy ion 
beams at the accelerator complex FAIR [11 . The con- 
ditions of high energy are beneficial from the theoretical 
point of view. First, under high-energy conditions the 
(p, 2p) interaction range is small compared to the size 
of the target nucleus and one can make use of the zero- 
range approximation when modeling the proton-nucleus 
interaction vertex. Second, for fast continuum nucleons 
with a sufficiently small de Broglie wavelength, the ef- 
fect of attenuation can be appropriately and accurately 
computed in the semi-classical eikonal approximation [12] 
even when adopting a Dirac treatment of the nucleons in 
the initial and final channel [T3]. From Fig. [l] it is clear 
that the proton de Broglie wavelength drops below 1 fm 
for nucleon kinetic energies larger than about 500 MeV. 
Several studies [14 have pointed towards the applicabil- 
ity of the eikonal method up to remarkably low nucleon 
kinetic energies of about 250 MeV. 

Just as for the A{e, e'p) reaction, the observables from 
A{p^ 2p) measurements are convoluted in that nuclear 
attenuation does not allow one to relate the measured 
ejected nucleon properties directly to the physics at the 
interaction point. With one proton subject to initial- 
state interactions (ISI), and two protons subject to FSI, 
the development of a reliable reaction theory is of the 
utmost importance for a quantitative analysis of the 
A{p,2p) data [H [I5l [16]. It is weh established that 
nuclear attenuation tends to make the detected signals 
less sensitive to the high-density regions of the target. 
Accordingly, one point of concern is whether or not one 
may learn something about the bulk properties of nuclei 
from A{p^2p). Or, in other words, carry the ejected nu- 
cleons information about the interior of the system or do 
they mostly originate from the peripheral regions of the 
nucleus? In this paper we wish to address this issue. We 
attempt to quantify what regions of the target nucleus 
can be probed in A{p^2p). For the sake of reference, 
we also add results for A{e^e'p). The latter serve as a 
benchmark, as with A(e, e'p) a lot has been learned about 
single-particle properties of stable nuclei during the last 
couple of decades. We focus on single- nucleon knockout 
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FIG. 1: The de Broglie wavelength (Xds) of a proton as a 
function of its kinetic energy. 



reactions with nucleon kinetic energies of several hun- 
dreds of MeV. At those energies the Glauber multiple 
scattering framework is both appropriate and accurate 
to model nuclear attenuation. 



In Sec. II A[ we outline the necessary formalism for 
quasifree single-nucleon knockout in a relativistic and 
cross-section factorized framework. In Sec. IIIBI a rel- 
ativized version of Glauber multiple scattering theory is 
introduced and it is pointed out how it can be corrected 
for the medium effects related to SRC. The results of the 
the numerical calculations are discussed in Sec. Illll and a 
summary is given in Sec. |IV[ 



II. THEORY 



Cross sections 



The theoretical calculations for the A{p^pN) and 
A(e, e'A^) reactions presented here are performed with 
a factorized form for the cross sections. They represent 
very useful zero-th order approximations and are often 
used as a starting point for the interpretation of the data. 
The factorization is largely based on the impulse approx- 
imation (lA) which develops out of a reaction picture 
in which {A — 1) nucleonic degrees-of- freedom are frozen 
during the interaction of the external probe and the tar- 
get. In the lA, the role of the spectator {A — 1) nucleonic 
degrees-of-freedom is restricted to distorting the waves of 
the impinging and ejected nucleons. In addition, the fac- 
torization allows for a more direct comparison between 
the (p, 2p) and the (e, e'p), as their differential cross sec- 
tions become proportional to a distorted momentum dis- 
tribution that is related to the probability of finding a 
nucleon with well-defined quantum numbers and a cer- 
tain momentum in the target nucleus. First, we derive 
a factorized expression for the A{e, e'p) and the A{p, 2p) 
cross section. 

Consider the A(e, e'N)A—l process and define the cor- 
responding kinematic variables of the impinging electron, 
the target nucleus, the scattered electron, the residual 
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nucleus, and the ejected nucleon as 

K^' (e, k) + K'X (Ea, Ua) (e', k') 



K>^_, (Ea-iSa-i) + K {En, kN) (1.) 



The fivefold differential cross section in the lab frame 
adopts the form [17 
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where P^a = | ~ 1 is the relative velocity and cor- 
responds to the appropriate average over initial states 
and sum over final states. We assume that all wave 
functions entering the computation of the reaction ampli- 



tude M 



(e,e'p) 
fi 



are normalized to unity. Throughout this 



work, the adopted normalization convention for the Dirac 
spinors is u (^k, rus^ u (^k, m^^ = 1. The recoil factor 



in Eq. ([2| reads 
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where the momentum transfer has been defined as 

q = k-k' . (4) 

For the kinematics of the A{p^ 2p)A — 1 we adopt the 
following conventions 

Pf (EpuPi) + K'X [EaSa) ^ K'l (EkiM) 

+ K!^(Ek2M)+K'X_^(EA-l,kA-l) ,(5) 

where refers to the impinging nucleon, and (i^f , ) 
to the pair of ejected nucleons. The fivefold differential 
equation for the 2p)A—l process can be straightfor- 
wardly derived from the A{e^ e'p) one (Eq. ([2|) by means 
of the following substitutions (mg Mat, e Epi^ k' 
ki^kN k2) and this results in 
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where the recoil factor reads 
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In this expression the momentum transfer is defined as 

q=pi-ki . (8) 

In the laboratory frame one obtains the following expres- 
sion for the relative velocity PpA = 



(p,2p) 
fi 



can be related 



The squared amplitude ^ M. 

to the free proton-proton cross section after making sev- 
eral assumptions. They are pointed out in great detail in 
Ref. [16] and include the neglect of the negative-energy 
components in the proton-proton scattering amplitude. 
Further, one assumes that the operator describing the 
collision of the impinging proton pi and the bound pro- 
ton pm resulting can be described in terms of the on-shell 
proton-proton scattering amplitude. Now, we sketch how 
one can arrive at a factorized expression for the ampli- 
tude TWj^'^^^ We denote the position coordinates of 
the impinging proton as fo, of the ejected protons as 
(n),^.), and of the nucleons in the target nucleus as 
(fi, . . . , f^). In the eikonal approach, the attenuation can 
be accounted for by means of a multiplicative factor to 
be applied to a plane- wave wave function. Accordingly, 
the relativistic distorted wave function of the impinging 
proton reads 



E^M 
2M 



1 

E+M 



where u{p^ s) is a four-component free-particle Dirac 
spinor. For the distorted wave functions 6? , and 

of the ejected protons, similar expressions hold. 



In the above expression, the ^d-body operator Spi is fully 
responsible for the effect of attenuation. In this work, we 
compute the effect of attenuation in a relativized Glauber 
mod el. M ore details will be provided in the next subsec- 
The Spi is a two-by-two matrix which acts on 
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tion 

the Pauli spinors. Here, we assume that Spi is a diagonal 
matrix which amounts to neglecting the spin-dependent 
attenuation mechanisms. Indeed, the central component 
of the nucleon-nucleon scattering amplitude accounts for 
the major impact of nuclear attenuation in high-energy 
proton- nucleus collisions [18 . Recent investigations [19] 
have clarified the role of the spin-dependent terms in the 
nucleon-nucleon scattering amplitude for the attenuation 
effects in D{e^ e'p)n. The central component was found 
to dominate the attenuation for most observables, and in 
particular for the differential cross sections at low miss- 
ing momenta. Accordingly, we deem that inclusion of the 
central component is sufficiently accurate for our cur- 
rent purposes, namely making a comparative study of 
the density dependence of the {e^e'p) and (p, 2p) reac- 
tion throughout the mass table. 

We now proceed with the derivation of a factorized ex- 
pression for the (e, e'p) and (p, 2p) cross section. In what 
follows we describe the wave function of the target nu- 
cleus by a normalized Slater determinant \a1a2 • • - o^a) 
where the ai refer to the quantum numbers of the oc- 
cupied single-particle states. With the distorted wave 
functions of Eq. ([9| we get after neglecting the negative- 
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energy projection term the following expression [16] 
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where Fpp is the pp scattering amplitude in momentum 
space. The quantum numbers {nnm) = ai determine 
the orbit of the struck nucleon which is described by a 
relativistic single-particle wave function ^n/^m(^)- The 
missing momentum p^ is determined by the difference 
between the asymptotic three- momentum of the ejected 
nucleon /c2 and the three-momentum transfer q = Pi — ki 



+ ^1 - Pi 



-kA- 



(11) 



In the absence of nuclear attenuation, the missing mo- 
mentum equals the momentum of the bound nucleon with 
quantum numbers {nnm) which collides with the proton 
beam. In the above equation ([lO| we have introduced an 
operator which accounts for the ISI/FSI [16^ 



1=2 

X Ski (r, r2, . . . , va) Sk2 (r, r2, . . . , va) {12) 



a multi-dimensional convolution over the squared wave 
functions of the spectator nucleons and a product of the 
scalar operators Sk for the impinging proton and two 
ejected protons 



With the aid of Eq. (10) one arrives at the following 
factorized form for the five-fold A{p^ 2p) differential cross 
section of Eq. (p|) 
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where f^SS") is the center-of-mass cross section for 

V ail /cm. 

pp scattering at an invariant energy W = ^/s. The Suk 
is the spectroscopic factor ( < Suk < 1) that is related 
to the occupancy of the orbit {nti) in the ground state of 
the target nucleus. In the above expression the distorted 



momentum distribution PnAPm) adopts the form [20] 

PuKiPm) = 22 ^^ (27r)3/^ ^^^"^^ ^)^RMSGA(0^n/.m(r) 

= Yl i^nnm {Pm)y ^uKm iPm) , 
s,m 

= j dr j d95 (r, 0) , 



where 



D{r,e) = J 



(14) 



d(j)r'^ sine u{pm, s)S^^^lj^{r)(j) 
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Working along similar lines one can derive a factorized 
form of the A{e^ e'p) differential cross section [21] 
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where cr^^ is the off-shell electron-proton cross section 
obtained from positive-energy projections. The distorted 
momentum distribution is defined as in Eq. ( 14 ) whereby 

ie.e'p) 
RMSGA 



•^RMSGA 



is replaced by S^-^<^qj^ that adopts the form [17] . 

i=A 



^RMSGa(^) = n / 1^"^ (^^)l^ '^^1 (^' ^2, • • • , Ta) . 
i=2 ^ 

(17) 



B. Relativistic multiple-scattering Glauber 
approximation 



The quantity 6{r^0)drd0 defined in Eq. (14) en- 



codes the contribution from an infinitesimal interval 
[r -\- dr^O -\- dO] to the cross section for a quasifree p-\-A 
p -\- p -\- A — 1 process that leaves the residual nucleus in 
a hole state determined by the quantum numbers (nhzm) 

[22]. The eikonal operator S^usGAi^) of Eq. ^ re- 
ceives contributions from the impinging (Spi) and the 

two ejected (Ski^Sk2) protons. In Refs. [17] and [23] 
we developed a relativized version of Glauber multiple- 
scattering theory. In this so-called RMSGA approach the 
eikonal phases are diagonal 2x2 matrices 



Spi(r{h,zyr2,...rA) = H [l - T (b - bj^ {zj - z) 

J=2 

= expix {r(b,z^ ,pi^ . (18) 
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For the profile functions TpN for proton-nucleon scatter- 
ing we adopt tiie standard Gaussian parametrization 



(1 - iCpN) 



exp - 



pN 



2/32 



pN 



wtiere •, PpN •, c^^^n have been determined from tiie 
database of proton-proton and proton-neutron cross sec- 
tions [17] . 

It can be expected ttiat tiie nucleon-nucleon interac- 
tions entering Eq. (18) by means of tiie profile function 
FpTv will be subject to medium modifications [24l [25] . 
Mechanisms like Pauli blocking often lead to an effec- 
tive reduction of the nucleon-nucleon cross sections in the 
medium. At higher energies, the effect of Pauli blocking 
is expected to become small. Another important source 
of medium effects in the treatment of ISI and FSI, are 
short-range correlations (SRC) [26H28]. The SRC are 
related to the finite size of the nucleons and the liquid 
properties of the nucleus. Indeed, the presence of a nu- 
cleon at some position r induces local fluctuations in the 
nuclear density. These local fluctuations, that go beyond 
the mean-field picture, can be included in the numerical 

evaluation of the operators 5^^MgQ^(r) and 5^^|Q^(r), 
where r is the interaction point of the impinging beam. 

With the Eq. ([l8|, the ^rm|ga(^) of Eq. ^ has a very 
intuitive interpretation: it represents the accumulated 
phase of one incoming and two outgoing waves which 
are subject to a medium with (A — 1) grey disks charac- 
terized by the profile function T^n and distributed over 
the medium by a density distribution of the mean-field 
type 111=2^ l^tti ifi)\^' The latter expression for the den- 
sity distribution can be corrected for SRC, by using the 
information that for a reaction to take place a nucleon 
should be present at the initial interaction point |29j. 
This can be achieved in the following way. First, the 



squared single-particle wave functions in the Eqs. (12) 



and (17) can be approximated by the one-body density of 



the target nucleus (r) (normalized as J drp^^ (r) = A) 



[1]/ 
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if I'/'o 



A 



(19) 



This substitution has a relatively minor impact on the 
computed effect of ISI/FSI in the RMSGA model [H]. 
Without any loss of generality the p^^ (r) can be replaced 
by the ratio of the two-body density p^^ (normalized as 

J dfi J dr2p^A(^^^^'^) ~ one-body den- 

sity: 



(20) 



where r is the coordinate of the probe-target interaction. 
For an uncorrelated two-body density one has: 



^A,uncorr. 



(^1,^2) 



A 



(21) 



and Eq. (20) becomes trivial. One can implement the 



effect of central SRC in the two-body density by adopting 
the following functional form [29l [30] 



,{ri,r2) = ^-^l{ri)p^^{ri)p^^{f2)l{r2)g{ri2) , 



(22) 

where g{r 12) is the so-called Jastrow correlation function 
[31 and 7(r) a function which guarantees the proper nor- 
malization of The 7(r) can be numerically obtained 
as the solution of an integral equation. With the above 
expression for the two-body density, the Eq. (20) be- 
comes 



P!l'(r2) 



l{r2)PA {r2)l{r)9{\r2-r\) 



pf{r2,r). (23) 



From the above derivations it follows that the computa- 
tion of ISI/FSI can be corrected for SRC by replacing 
\(t>a. (ri)l^ with pf{ri,r)/A in the Eqs. ^ and 
The presence of a nucleon at the interaction point r in- 
duces local fluctuations in the probability distributions of 
the remaining A — 1 nucleons. Within a radius r < r^^ 
where tn is the radius of a nucleon, of the initial interac- 
tion point, there will be a reduced probability to find a 
nucleon to scatter from, whereas for r ^ 2rN there will be 
an enhanced probability. The sole input required to de- 
termine p^(r2,r) from p{^^(r2) is the Jastrow correlation 
function g{r 12). As this function is related to the short- 
range dynamics of nuclei it is considered "universal" . We 
use a g{r 12) that has a hard core extending over 0.8 fm 
and a second bump with a peak at ri2 ~ 1.3 fm. With 
this choice for the g{r 12) we obtained a fair description 
of the SRC contribution to the exclusive {e^e'pp) cross 
sections from [32 and ^^O [33 . In Fig.[2]we display 
for ^^C and ^^^Pb the effective density for a nucleon that 
is hit by an external probe in the center of the target nu- 
cleus (z=0 fm) and in a more peripheral location along 
the z-axis at 2 fm from the center {z = 2 fm) of the tar- 
get nucleus. The SRC induce a hole in the density at the 
position of the probe-target interaction point, and some 
enhancement at distances ~ 1.5 fm further away. The 
proposed method for implementing the effect of SRC in 
Glauber calculations leads to effective densities (Fig. [2| 
that are qualitatively very similar to those produced in 
ab-initio calculations (see for example Fig. 1 of Ref. [34^ ) . 
Therefore, we consider the proposed method to account 
for SRC in the ISI/FSI calculation, as realistic and effi- 
cient. 



III. RESULTS 

We now present the results of the numerical calcula- 
tions for S{r^O). The function S{r^O) depends on the 
quantum numbers of the bound nucleon that collides with 
the proton or electron beam and the kinematics of the re- 
action. The considered A(e, e'p) and A(p, 2p) differential 
cross sections involve five kinematic variables. The pre- 
sented investigations aim at determining global trends 
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C MF 

C MF+SRC ^ = 
C MF+SRC z = 2 
Pb MF 

Pb MF + SRC 2=0 
Pb MF + SRC 2=2 




*5rmsga = 1 in Eqs. (12) and (17). Under those con- 



0.00^ 



FIG. 2: (Color online) Comparison of the mean- field (MF) 
density p^^{r) and the SRC corrected effective density (de- 
noted as MF+SRC) p^/^ {r, r = (0, 0, z)) for the target nuclei 
^^C and 2°^Pb. The dashed (dotted) hnes are for ^ = 
(^ = 2). 





FIG. 3: Schematic representation of the A[e^ep) reaction 
in parallel kinematics and the A{p^ 2p) reaction in coplanar 
and symmetric kinematics. The quantity 6{r, 0)drdO provides 
the contribution from the interval [r + c/r, ^ + dO] to the cross 
section. 



and choices with regard to the kinematics are in order. 
For the A(e, e'p) reaction we restrict ourselves to so-called 
parallel kinematics: the final nucleon is detected along 
the direction of the momentum transfer q. For A(p, 2p) 
we consider coplanar and symmetric kinematics: the two 
ejected nucleons have equal kinetic energies (| hi | = | ^2 |) 
and escape with equal opening angle but on opposite 
sides relative to the direction of the momentum trans- 
fer. The considered kinematics is illustrated in Fig. |3] 
We stress that in many respects there are strong analo- 
gies between the selected kinematics so that meaningful 
comparisons between (e, e'p) and (p, 2p) can be made. 
The effect of ISI and FSI can be nullified by setting 



ditions one adopts the impulse approximation and all 
impinging and ejected nucleons are described by plane 
waves. Accordingly, we will refer to the corresponding 
reaction model as the relativistic plane- wave impulse ap- 
proximation (RPWIA). The difference between the RP- 
WIA and RMSGA predictions for J(r, 0) can be exclu- 
sively attributed to nuclear attenuation. In the calcu- 
lations including attenuation, we discriminate between 
the "RMSGA" and the "RMSGA+SRC" model variants. 
The latter includes the effect of SRC in the modeling of 
the ISI/FSI. 

In what follows we will display (e, e'p) and (p, 2p) re- 
sults for the target nuclei "^He, ^^C, ^^Ye and ^ospb and 
use the relativistic cfuj model to determine their single- 
particle wave functions. We define the 2:-axis along the 
q and the xz-plane as the reaction plane. The ejected 
nucleons are detected in the forward direction. We will 
start our discussions with the results for the medium- 
heavy target nucleus ^^Fe. 

In Fig. [4] we display the function (5(r, 6) for proton 
knockout from the 1^1/2 , 1^3/2 7 l<^3/2 7 and I/7/2 orbits 
from ^^Fe. We compare the (p, 2p) with the (e, e'p) result 
for an energy transfer of 1.5 GeV and conditions probing 
the maximum of the undistorted momentum distribution 



{Pm) 



E 



dr 



(2^)3/2 



u{pm,s)(j)nKm{r) 



(24) 



The X-axis is a symmetry axis for 5(r, 0) in the RPWIA 
reaction picture. Indeed, in the absence of nuclear atten- 
uation, the upper (0° < 6> < 90°) and lower hemisphere 
(90° < 9 < 180°) of the target nucleus equally contribute 
to the measured signals in the detectors and the (5(r, 6) 
becomes equal for {e^e'p) and (p, 2p). In RPWIA the 
(5(r, 0) reflects the symmetry imposed by the quantum 
numbers {nhz). It is obvious that with increasing orbital 
angular momentum / the detected signals are increasingly 
stemming from the peripheral areas of the target nucleus. 
The ISI and FSI have the strongest impact at the highest 
nuclear densities. This reflects itself in the fact that the 
largest values of S{r^ 0) in the RMSGA model are shifted 
to larger values of r in comparison with what one finds 
in RPWIA. As can be appreciated from Fig. |4] this shift 
is strongest for the deep- lying lsi/2 orbit and less pro- 
nounced for the valence I/7/2 orbit. Further, attenuation 
breaks the symmetry between the contribution from the 
upper and lower hemisphere as the first one is positioned 
closer to the detector (s). The degree of asymmetry be- 
tween the contribution to the cross section from (r, 0) 
and (r, 180^ — ^) is a measure for the impact of attenua- 
tion. Another indicator is the ratio of the magnitude of 
RMSGA to the RPWIA prediction for S{r,0). Loosely 
speaking one could interpret this ratio as the fraction 
of the available signal in the target nucleus at some po- 
sition (r, ^) that can withstand the attenuating nuclear 
medium and makes it to the detectors. Obviously, the 
induced angular asymmetry, radial shift and overall re- 
duction occur for the S{r^O) in both (e^e'p) and {p,2p). 
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FIG. 4: (Color online) The function 6{r, 0) for knockout from various shells in the ^^Fe target nucleus and ejected nucleon kinetic 
energies of 1.5 GeV. The (e, e'p) results are for parallel kinematics. For (p, 2p) coplanar and symmetric kinematics is considered. 
The magnitude of the momentum transfer |^ is adjusted so as to probe the maximum of the momentum distribution. This 
corresponds with prn = MeV for the lsi/2, Pm = 105 MeV for the lp3/2, Pm = 145 MeV for the lc/3/2, and pm = 180 MeV 
for the 1/7/2- For the sake of reference, the measured proton root-mean-square radius in ^^Fe is Vrms ~ 3.75 fm ^35^ . 



All three effects, however, are far more pronounced for 
the (p, 2p) than for the corresponding 5{r^0) in {e^e'p). 
Further, one observes that the combined effect of ISI/FSI 
gradually diminishes as one moves from the deep-lying to 
the valence single-particle orbits. 

The preceding discussion concerned one particular pro- 
ton kinetic energy. Next, we report on the study of the 
energy dependence of the impact of nuclear attenuation. 
As a representative example, in Fig.jsjwe show for (e, e'p) 
and (p, 2p) the radial reaction probability distribution 
5{r) = J d05 (r, 0) for proton knockout from ^He at ki- 
netic energies ranging from 0.25 to 1.8 GeV. One observes 
a very soft energy dependence in the radial dependence 
and the magnitude of the 5{r). Therefore, we deem that 
the results for the average densities and radii that will be 
presented below and that are obtained at a specific ki- 
netic energy of 1.5 GeV can be considered as representa- 
tive for nucleon knockout reactions with T/v ^ 250 MeV. 



In Figs. [5)[6)and[7| we have added the function r'^p\^ (r). 
The densities are computed with the basis of relativistic 
single-particle wave functions which is also used in the 
reaction dynamics calculations. The densities are plot- 
ted not to scale and help in evaluating to what extent the 
knockout process succeeds in probing the bulk regions of 
the target nucleus. 

We wish to investigate the target-mass dependence of 
5{r^0). The results of Fig. [4] pointed towards a strong 
orbit dependence in the average radius and density that 
can be probed. We consider knockout from the deep- 
lying 1^1/2 (Figj6| and one of the weakly-bound orbits 
(Figj7| for the four representative target nuclei consid- 
ered in this work. For ^^^Pb we opted for the l/iii/2 
valence orbit as it is representative for orbits with large 
orbital momentum. The I/111/2 ^ hole state in ^^'^Tl is 
located at an excitation energy of 1.35 MeV. Results of 
the RPWIA and RMSGA calculations are contained in 
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FIG. 5: (Color online) The energy dependence of the radial 
reaction probability densities 5{r) — J d06{r,0). We consider 
proton knockout from ^He at ejected proton kinetic energies 
of 250, 400, 600, 900, 1200, 1500, and 1800 MeV. The (e, ep) 
((p, 2p)) RMSGA results are for parallel (coplanar- symmetric) 
kinematics. For the sake of reference the r^p^^^ (r) for ^He is 
shown (not to scale!). 
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FIG. 6: (Color online) The target mass dependence of the 
radial reaction probability density S(r). We consider proton 
knockout from the deep- lying lsi/2 level and ejected proton 
kinetic energies of 1500 MeV. The magnitude of the momen- 
tum transfer |^ is adjusted so as to probe the maximum of 
the momentum distribution {pm = MeV). For the sake of 
reference we added pa(t) for the various target nuclei (not 
to scale!). 
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FIG. 7: (Color online) As in Fig. [6]but now for knockout from 
one of the valence orbits. 



the Figs, [g] and [tI We compare the (5(r) for the (p, 2p) re- 
action with the (e, e'p) one. The RPWIA result for 5(r) 
is the reference figure of merit for knockout from a par- 
ticular orbit. The mass dependence can be best appre- 
ciated from the 1^1/2 knockout results from Fig. |6] The 
most spectacular observation is the enormous decrease in 
the absolute value of the RMSGA (p, 2p) radial reaction 
transition density for medium-heavy and heavy nuclei. 
For the 1^1/2 orbit, attenuation sheds about 90% of the 
(p, 2p) RPWIA signal in ^^C, for ^^^Pb this becomes more 
than 99%. It is clear that for any meaningful extraction 
of spectroscopic factors from (p, 2p) data the availability 
of a highly reliable reaction model is of the utmost im- 
portance. For knockout from the valence orbits (Fig. [t]) 
the major fraction of the signal stems from the tail of 
the density distribution of the target nucleus. For the 
light nuclei, the (e, e'p) reaction performs only slightly 
better than (p, 2p) when it comes to probing the bulk 
regions. For the heavy nuclei, the growing role of attenu- 
ation makes the reaction probability density increasingly 
peripheral, and it is obvious that this primarily affects 
the (p, 2p). 

In order to quantify the average density and radius that 
can be probed in the quasifree single-nucleon knockout 
reactions, we introduce the quantities [22l [36] 



(r) 



/ drddp^^} (r) 8 (r, Q) 

/ drdBb (r, B) 
/ drddrS (r, 6») 
/ drde8 (r, Q) ' 



(25) 
(26) 



where (r) is the density of the target nucleus. From 
the foregoing discussions we infer that for a given tar- 
get nucleus A and orbit, the energy dependence of the 
effective densities (p) and effective radius (r) is rather 
smooth. For that reason we stick with Tp=1.5 GeV and 
compile in Figs, [s] and [9] the (p) and (r) for the (e, e'p) 
and (p, 2p) results contained in Figs, [g] and [7| The RP- 
WIA prediction for the average radius (r) increases with 
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FIG. 8: (Color online) The target-mass dependence of the average density (p) and the average radius (r) that can be probed 
in single-nucleon knockout from the deep-lying Is shell. The kinematic conditions are those from Fig. [6] The dotted line is 
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FIG. 9: (Color online) The target-mass dependence of the average density {p) and the average radius (r) which can be probed 
in single-nucleon knockout from the valence shell. The kinematic conditions are those from Fig.[7| The dotted line is 1.2A^^^. 



A at a slightly softer rate than A^^^. Note that even 
in the idealized attenuation-free world described by RW- 
PIA, there is a strong mass dependence in the average 
density that can be probed in a quasifree nucleon knock- 
out process. This strong mass dependence was exploited 
in Ref. [22] in order to probe the medium dependence of 



the A^A^ scattering amplitude. In comparing extracted 
information from (e, e'p) and (p, 2p) reactions, like spec- 
troscopic factors, through the mass table, it is often not 
realized that the reactions probe increasingly the target's 
surface region with growing A. For proton knockout from 
^He, RPWIA predicts an average density {p) = O.SSpo, 
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FIG. 10: (Color online) The RPWIA and RMSGA distorted 
momentum distribution for knockout from the deep-lying 
lsi/2 orbit. 
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FIG. 11: (Color online) The RPWIA and RMSGA distorted 
momentum distribution for knockout from the valence shell. 



with the nuclear saturation density po = 0.16 fm~^. In 
^^^Pb and knockout from the I/111/2 valence shell we find 
(p) = O.OSpo- All this before correcting for FSI/ISI which 
further reduces the (p). 

For light nuclei the RPWIA and RMSGA prediction 
for (r) are close. For heavier nuclei a different story 
emerges. For knockout from the lsi/2 orbit in ^^^Pb 
for example, attenuation makes the (r) to grow quite 
dramatically. Indeed, For {p,2p) ((e,e'p)) the RMSGA 
prediction is ~ 2.5 fm ( ~ 0.9 fm) larger than the RP- 
WIA value of 4.71 fm. For the valence states (Fig. [9| 
the increase in (r) is not larger than 0.4 fm in (e, e'p). 
In (p, 2p) the gain in (r) is larger, but even in Pb it is 
smaller than 1 fm. As mentioned before, the (p) is de- 
creasing monotonically from the lightest to the heavier 
nuclei. Attenuation makes that in reality smaller (p) will 
be probed. The effect varies from a loss of couple of per- 
cent to a loss in (p) of more than 50%. In ^^C, the (p, 2p) 
reaction from the si/ 2 orbit can effectively probe higher 
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FIG. 12: The predicted average density which can be probed 
in a quasifree nucleon knockout reaction. The results are ob- 
tained in the RMSGA framework and include the effect of 
SRC. 



densities ((p) = 0.39po) than the {e^e'p) reaction from 
the valence P3/2 orbit ((p) = 0.25po)- For the -^^C(p, 2p) 
reaction with knockout from the 5l/2 orbit the RMSGA 
prediction of {p) = 0.39po is comparable to the DWIA 
result (p) = 0.34po of Ref. [36 which is obtained for 1 
GeV incoming protons. This illustrates the robustness of 
the results of this work. 

The results for (p) and (r) in Figs, [s] and |9] allow one 
to estimate the role of SRC in the modeling of the at- 
tenuation. In line with the observations of Refs. and 
[30] the SRC make the nucleus somewhat more trans- 
parent. This reflects itself in the fact that after includ- 
ing the SRC the RMSGA predictions for (p) and (r) are 
shifting towards the RPWIA values. We wish to empha- 
size that the results of Figs. [S] and |9] refer to kinematics 
corresponding with the maximum of the single-particle 
momentum distributions. On the basis of the densities 
shown in Fig. |2] one may be tempted to expect rather 
spectacular effects from SRC in the ISI/FSI. The effect 
of SRC on the (p) and (r) is rather moderate due to the 
fact that the ISI/FSI are long-ranged in the longitudi- 
nal direction and short-ranged in the transverse direction 

[26;. 

Up to this point we have evaluated quasifree processes 
for which the kinematics is tuned to probe the maxi- 
mum of the momentum distribution. Now, we turn to a 
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study of the distorted momentum distribution p^niPm) as 
a function of the missing momentum pm- We stick with 
parahel kinematics for (e, e'p) and coplanar and symmet- 
ric kinematics in (p, 2p) . The variation in pm is achieved 
by varying the kinetic energies of the eject iles at a fixed 
value of the momentum transfer q. We study knock- 
out from the deep- lying lsi/2 level for the various nuclei 
in order to get some feeling about the mass dependence 
of the attenuation. As one can appreciate from Fig. [lO] 
with growing A the RPWIA and RMSGA predictions for 
PnK.{Pm) increasingly diverge. A similar remark applies 
to the pnK{Pm) for {e^e'p) and (p, 2p). For the valence 
states (Fig. [IT]) one observes similar trends, though less 
pronounced. This can be easily understood by consid- 
ering that they probe the peripheral areas of the target 
nucleus. 

At large missing momenta Pm the effect of the SRC cor- 
rections on the distorted momentum distributions can be 
extremely large [26 . Here, we concentrate on low missing 
momenta and observe that SRC make the nucleus more 
transparent for the emission of nucleons. This observa- 
tions complies with the conclusions of Refs. [26] and [30] . 
The SRC do not dramatically affect the Pm dependence 
of the distorted momentum distributions at low Pm but 
have an effect on its magnitude. This makes it of the 
utmost importance to use a correlated Glauber approach 
in order to extract precise information about the spec- 
troscopic factors SriK- 

IV. CONCLUSIONS 

In summary, we have used a relativistic framework to 
make a comparative and consistent study of the effective 
nuclear densities which can be probed in quasifree (p, 2p) 
and (e, e^p) reactions at high energies. We use relativistic 
single-particle wave functions from the auj model and a 
relativistic extension of Glauber multiple-scattering the- 
ory. Adopting the impulse and factorization factoriza- 
tion, both the (e, e^p) and the (p, 2p) cross sections are 
proportional to the distorted momentum distribution 
which include the effect of nuclear attenuation for the im- 
pinging and ejected nucleons. The reflects the effec- 
tive momentum density for a bound nucleon with quan- 
tum numbers nHi that is accessible for a certain reaction. 
We use this quantity as a tool to make a quantitative as- 
sessment of the role of nuclear attenuation for both types 



of reactions. Obviously, nuclear attenuation is a compli- 
cating factor in the extraction of nuclear-structure infor- 
mation from the measured single-nucleon knockout sig- 
nals. Strong attenuation has a severe geometric influence 
as it makes the detected signal to carry little informa- 
tion about the nuclear interior. An interesting question, 
therefore, is how efficient (p, 2p) reactions are in probing 
the nuclear interior in comparison with the time-honored 
{e.e'p). 

We have presented numerical results for knockout from 
one of the weakly-bound levels and knockout from the 
deep-lying Is level from "^He, ^^C, ^^Fe, and ^ospb. Jt 
emerges that to a remarkable extent the effect of the 
nuclear attenuation on the angular cross sections is in- 
dependent of the energy of the initial and final protons. 
Accordingly, we have focused on one energy which we 
consider representative. The results for the average den- 
sities (p) which can be probed in quasifree single-nucleon 
knockout are collected and shown in Fig. [12] The (p)'s 
are very sensitive to the quantum numbers of the bound 
proton which collides with the probe. Further, there is 
strong target-mass dependence. Whereas in the lightest 
nuclei one can probe densities comparable to nuclear sat- 
uration density po, for a mid- heavy nucleus like Fe this is 
of the order of 10% of po and even smaller average den- 
sities are probed for a heavy nucleus like Pb. For light 
nuclei, the (p, 2p) and (e, e'p) reactions are comparably 
efficient in probing the nucleus' inner regions. In ^^^Pb, 
the effect of attenuation is very dramatic for (p, 2p) and 
the average density which can actually be probed can be 
half of the (e, e'p) one. 

Our model for the ISI and FSI implements the effect 
of short-range correlations. In line with previous stud- 
ies we find that the SRC tend to reduce the influence 
of attenuation. The SRC moderately affect the average 
density and radius that can be probed. We do find, how- 
ever, that after correcting for SRC the cross sections can 
become substantially larger. The SRC corrections are 
particularly relevant for the heavier target nuclei and the 
channels which probe the inner-bound nucleons. As spec- 
troscopic factors are typically obtained from the ratio of 
the measured and the computed distorted momentum 
distribution, the SRC should become an essential ingre- 
dient of any model for nuclear attenuation. 

This work was supported by the Fund for Scientific 
Research Flanders. 
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